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A  system  for  ethanol  steam  reforming  and  purification  of  carbon  monoxide  (CO)  designed  to  feed  a  PEM 
fuel  cell  has  been  modelled.  From  the  model,  we  study  the  sensitivity  and  controllability  emphasizing 
the  study  of  the  influence  of  the  temperature  on  the  output  variables  of  interest.  The  results  of  the  study 
of  controllability  are  used  for  the  identification  of  the  best  control  structures. 
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1.  Introduction 

The  gradual  reduction  of  fossil  fuel  reserves  and  the  envi¬ 
ronmental  pollution  problems  associated  with  their  combustion 
have  turned  the  attention  of  researchers  to  the  search  for  alterna¬ 
tive  energy  carriers.  The  energy  vector  that  is  currently  receiving 
more  attention  is  hydrogen,  considered  as  a  possible  candidate 
to  partially  replace  oil  as  fuel  in  mobile  applications.  However, 
hydrogen  can  only  be  considered  a  clean  alternative  if  it  comes 
from  clean  energy  sources.  Ethanol  can  be  considered  a  source  of 
H2  neutral  in  relation  to  CO2  emissions.  This  has  motivated  many 
investigations  focusing  on  the  design  of  reactors  for  the  produc¬ 
tion  of  H2,  but  less  attention  has  been  given  to  the  control  of 
the  developed  systems,  that  includes  tasks  such  as  the  selection 
of  the  control  structures  and  the  design  and  tuning  of  the  con¬ 
trollers.  There  exist  different  tools  for  the  task  of  selecting  the 
controlled,  manipulated  and  measured  variables  and  link  these 
variables  to  make  control  loops  [2].  The  main  mathematical  meth¬ 
ods  for  designing  control  structures  are  based  on  the  relative  gains 
[2]  and  singular  value  decomposition  analysis  [1,2].  In  this  work 
we  select  the  most  appropriate  control  structures  for  an  ethanol 
reformer. 
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Reforming  process 

-  Ethanol  dehydrogenation  to  acetaldehyde  and  hydrogen  (EtOHD) 

-  Acetaldehyde  reforming  (AcR) 

-  Water  gas  shift  reaction  (WGS) 

Purification  unit 

-  Condenser  (CON) 

-  Preferential  oxidation  of  CO  (CO-PrOx) 

2.  Description  of  the  ethanol  steam  reforming  process  and 
the  dynamic  model 

In  previous  publications,  we  have  described  a  non-linear 
dynamic  model  for  a  low-temperature  ethanol  steam  reformer 
based  on  a  cobalt  catalyst  [5].  It  consists  in  a  tubular  reactor  loaded 
with  catalytic  monoliths.  The  reforming  process  is  divided  in  three 
separated  stages  (EtOHD,  AcR  and  WGS).  Experimental  data  for  the 
kinetics  of  the  two  first  stages  was  given  in  [3].  In  this  work,  we 
have  considered  in  addition  a  purification  unit,  which  is  required 
to  directly  provide  hydrogen  from  the  reactor  to  a  fuel  cell. 

2. 1 .  Description  of  the  ethanol  steam  reforming  process 

2.1.1.  Reforming  unit 

In  the  first  stage,  ethanol  reacts  over  Sn02  to  produce  acetalde¬ 
hyde  and  hydrogen  according  to  the  reaction: 

c2h5oh  -*  C2H40  +  H2  (1) 
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Nomenclature 

C  concentration  (mol  mR-3 ) 

cp  specific  heat  (J  kg-1 I<-1 ) 

dt  monolith  external  diameter  (m) 

F  molar  flowrate  (mol  s-1 ) 

k  reaction  rate  constant  (mLgcat_1  h-1 ) 

I<  equilibrium  constant  for  homogeneous  reaction 

L  reactor  length 

t  time  (min) 

T  temperature  (I<) 

Tf  furnace  temperature  (K) 

U  overall  heat-transfer  coefficient  (J  s-1  m-2  k-1 ) 

y  mole  fraction 

p  pressure  (atm) 

v  superficial  velocity  (ms-1) 

x  conversion 

r  velocity  rate  (mol  mR-3  s-1 ) 

?7h2  hydrogen  yield  (%) 

OP  operating  point 

Subscripts 

e  equilibrium 

in  reactor  input 

out  reactor  output 

gas  gas 

i  reaction  number,  i  =  (ethanol  decompo¬ 

sition,  acetaldehyde  reforming,  water-gas  shift, 
respectively) 

j  component  number,  j  =  l,...,6  (C2H5OH,  H20, 

C2H40,  H2,  CO,  C02,  respectively) 

C2H5OH  relative  to  ethanol 
H20  relative  to  water 
C2H40  relative  to  acetaldehyde 
H2  relative  to  hydrogen 

CO  relative  to  carbon  monoxide 

02  relative  to  oxygen 

s  solid 

51  stage  1 

52  stage  2 

53  stage  3 

Greek  letters 
A  increment 

Vy  stoichiometric  coefficient  of  component  j  in  reaction 

i 

p  density  (kg  m-3) 


The  order  of  the  reaction  was  established  in  [3]  by  changing  the 
ethanol  load  at  a  series  of  temperatures  between  573  and  673  K.  A 
first-order  reaction  was  identified  and  a  kinetic  law  provided. 

The  acetaldehyde-steam  mixture  is  transformed  over 
Co(Fe)/ZnO  catalyst  in  the  second  stage  into  a  mixture  of  hydrogen, 
carbon  monoxide,  and  carbon  dioxide,  according  to  the  reactions: 


c2h4o  +  h2o 

3H2  +  2CO 

(2) 

c2h4o  +  h2o 

5H2  +  2C02 

(3) 

CO  +  h2o  -> 

co2  +  h2 

(4) 

For  the  kinetics  of  reactions  (2)  and  (3)  a  dependency  on  only 
acetaldehyde  concentration  was  encountered  [3].  The  kinetics  of 
the  reaction  (4)  is  based  on  the  partial  pressures  of  all  the  compo¬ 
nents  involved  in  the  reaction.  The  velocity  rates  of  these  reactions 


are: 


rC2H40  =  fcl,2  eEa-c2H4o/RTCc2H40  (5) 

fCO,S2  =  /<3  e£a-c°/RTPcoPH20  (l  -  -j? - :CT2  I  (6) 

\  Keq3  PcoPh20  J 

The  values  of  the  reaction  constants  (/<i  — /<3 )  are 
1 .38  x  1 04  mLgcat-1  h-1 ,  1.52  x  104  mLgcat_1  h-1  and  9.2  xlO2 
mol  m-3  s-1  atm-2,  respectively.  The  activation  energies 
(£a,c2H4o,  Fa  co )  for  these  equations  are  9.84  x  1 04  J  mol-1  and 
1.2  x  105Jmol-1,  respectively  [3]. 

In  the  third  stage,  a  Fe203-Cr203  catalyst  has  been  used  to 
further  carry  out  the  WGS  for  decreasing  the  CO  content  in  the  refor¬ 
mate.  In  this  case,  the  experimental  data  are  taken  from  [4].  This 
work  has  been  chosen  because  it  is  based  on  operating  conditions 
(rank  of  temperatures  and  pressure)  similar  to  the  conditions  of  our 
system.  In  this  case,  the  Langmuir-Hinshelwood  (L— H)  model  (7) 
describes  adequately  the  reaction  behaviour  over  the  temperature 
and  concentration  ranges  investigated  [4]. 


rCO,S3  =  kcO,S3  fCcO  f^H20 


Pc  o  Ph2o  (Pco2  Ph2  )/ Ke 


1  +  Kc  o  Pc  o  +  Kh2o  Ph2o  +  ^h2oPh2o 


(7) 


2.2.  Purification  unit 


The  purification  unit  has  been  specifically  included  in  this  work 
and  comprises  a  condenser  and  a  CO  preferential  oxidation  (CO- 
PrOx)  reactor  based  on  a  Pt  catalyst  [6,7].  The  condenser  is  used  to 
cool  hot  vapors  at  the  output  of  the  third  stage  and  separate  hydro¬ 
gen  and  carbon  oxides  from  excess  water  and  unreacted  ethanol 
and  acetaldehyde,  whereas  the  purpose  of  the  PrOx  reactor  is  to 
reduce  the  concentration  of  CO  at  the  inlet  of  the  fuel  cell  stack.  In 
the  PrOx  reactor,  air  is  mixed  with  the  feed  stream  to  preferentially 
oxidize  CO  (8)  while  minimizing  the  consumption  of  FI2  (9): 

CO+lo2^C02  (8) 

H2  +  1o2H20  (9) 

As  indicated  in  [7],  the  reaction  rates  in  the  CO-PrOx  stage  are 
described  by: 


rco  =  2/<Q2xo2 

Y^o^  +  fcco  y/*o7 
\/Xh2  2/<o2*o2 

f  h  = - 

2  fee OV^CO  +  fee  OV^CO 

2fo2  =  rco  +  fH2 
where 


fco2 


=  8.9  x  106  e10700/(RTs) 


Kco=  13.7 


(10) 

(11) 

(12) 

(13) 

(14) 


To  simulate  the  condenser,  we  have  solved  a  system  of  non¬ 
linear  equations  with  MATLAB™. 


2.3.  Dynamic  model 

To  obtain  the  dynamic  model  of  the  system,  pressure  is  con¬ 
sidered  constant  at  one  atmosphere  and  the  volumetric  velocity 
variation  along  the  reactor  is  taken  into  account. 
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Fig.  1.  Reforming  and  purification  process  of  an  ethanol-water  mixture. 


In  the  purification  unit,  isothermal  conditions  are  assumed. 
Radial  variations  are  neglected  and  only  the  axial  profiles  are  con¬ 
sidered. 

The  mathematical  model  is  based  on  the  mass  balance  and  the 
energy  balance,  and  it  is  assumed  that  the  reactor  is  a  pseudoho- 
mogeneous  system  by  considering  that  the  gas  temperature  and 
the  solid  temperature  are  the  same  [8]. 

Mass  balance  (for  component  j): 


8C  dvC 


if  T r» 


Energy  balance: 


(15) 


dr  dr 

Pgas  CPgas  +  Ps  CPs  -fo  +  Pg  ^gas1^ 


(16) 


The  energy  balance  parameters  [9]  are  17=4000  [J  s-1  K-1  m-3], 
ps  =  500  [kg  nrr3  ]  and  the  density  of  the  gas  has  been  obtained  from 
the  following  equation: 


£  _  PMmi 

dgas  -  - 

1X1  gas,  si 


(17) 


To  simulate  the  dynamic  behaviour,  initial  conditions  and 
boundary  conditions  are  set  as: 

Initial  conditions 


Q(0, x)  =  Qo(x)  with  xe[l,  L],j  =  1, 2, ... 6 
T(0,x)  =  T0(x)  with  xe[l,L] 


Boundary  conditions 

Cj(t,  0)  =  Cje(t)  with  t  >  0 
T(t,  0)  =  Te{t )  with  t  >  0 


the  goal  is  to  perform  a  sensitivity  and  controllability  analysis  of  the 
reforming  process  emphasizing  the  influence  of  the  temperature 
variables.  We  have  considered  as  inputs  the  flowrates  of  ethanol 
and  water  at  the  reactor  entrance  (Fc2h5oh,  FH2o),  the  temperature 
of  the  entering  mixture  (Tg  in),  and  the  temperatures  of  the  furnaces 
of  the  three  reforming  stages  (Tfsl,  Tfts2»  Tfs3).  As  outputs  (control 
objectives),  we  have  selected  the  flowrates  of  H2  and  CO  (Fh2  and 
Fco)  at  the  output  of  the  reformer. 

Fig.  1  shows  the  multiple  inputs  (manipulated  variables)  and 
multiple  outputs  (control  variables)  of  the  MIMO  system. 

The  linear  model  in  state  space  has  the  form: 

x  =  Ax  +  Buu  ,  m 

z  =  Cx  +  Duu  1  J 


where  the  state  x  contains  the  concentrations  of  all  components  in 
the  different  volumes  what  the  reactor  has  been  divided.  The  input 
vector  u  contains  the  manipulated  variables  and  the  output  vector 
z  contains  the  controlled  variables. 


^  =  [Fh2,  Fco] 

u  =  [Fc2H5OH  Fh20  Tg.in  TV, SI  Tf,S2  ^f,S3  ]T 


(21) 


2.5.  Linear  and  non-linear  models  comparison 


In  this  study,  the  nominal  steady-state  is  selected  at  a  H2  yield 
value  of  ?7h2  =  80  %  and  a  molar  fraction  of  carbon  monoxide  of 
yCo  =  0.83%  (see  Table  3),  which  is  suitable  for  entering  the  CO-PrOx 
reactor.  The  non-linear  model  is  evaluated  at  four  different  operat¬ 
ing  points  that  correspond  to  the  following  variations  with  respect 
to  the  nominal  operating  point: 


The  numerical  solution  of  the  partial  differential  equations  (15) 
and  (16)  was  accomplished  by  its  transformation  into  an  ODE- 
system  by  discretization  of  the  spatial  derivative.  To  this  end, 
backward  finite  differences  have  been  selected  (first-order,  15  dis¬ 
cretization  points  for  the  reported  simulations)  for  the  different 
stages  of  the  reforming  unit.  The  resulting  285  ODE  equations 
were  solved  by  an  algorithm  implemented  in  MATLAB™  (ODE45 
Normand-Prince).  Additional  details  regarding  the  mathematical 
model  can  be  found  in  [5]. 

2.4.  The  linear  model 

A  linearized  model  was  obtained  from  the  non-linear  model  of 
the  system  using  SIMULINK®  linearization  tools,  as  reported  in  [3]. 
The  state  space  representation  of  the  system  has  a  very  large  dimen¬ 
sion  and  therefore,  a  model  reduction  strategy  has  been  applied  to 
transform  the  original  model  into  a  simplified  form  that  has  lower 
order  and  preserves  the  dynamic  characteristics  of  the  original 
high-order  system  [5].  The  linear  model  is  the  base  of  the  control¬ 
lability  study  developed  in  this  work.  In  addition,  the  comparison 
between  the  non-linear  and  the  linear  models  will  provide  informa¬ 
tion  about  the  suitability  of  selecting  linear  controllers.  In  this  work, 


•  OPi  (AFc2h5oh)‘  10%  ethanol  input  increase  while  the  other  five 
inputs  are  kept  at  the  nominal  values. 

•  OP2(AFh2o)-  10%  water  input  increase  while  the  other  five  inputs 
are  kept  at  the  nominal  values. 

•  OP3  ( ATg  in):  10%  gas  temperature  input  increase  while  the  other 
five  inputs  are  kept  at  the  nominal  values. 

•  OP4  (ATfS2):  5%  furnace  temperature  of  stage  2  input  increase 
while  the  other  five  inputs  are  kept  at  the  nominal  values. 

The  response  of  Fh2  following  step  changes  in  the  different 
inputs  are  shown  in  Fig.  2,  where  both,  the  linear  and  the  non¬ 
linear  model  results  have  been  plotted.  It  can  be  observed  that  the 
linear  model  curves  have  minimal  differences  with  respect  to  the 
non-linear  model  curves  in  the  case  of  changes  in  FC2h5oh>  Fh2o> 
and  Tg  in.  On  the  contrary,  the  difference  is  notable  when  there  is  a 
change  in  Tf  S2.  This  non-linearity  can  complicate  the  control  of  the 
system  if  TfS2  is  used  as  a  manipulated  variable  to  control  Fh2. 

Fig.  3  shows  the  Fc o  output  response.  The  differences  between 
the  linear  and  non-linear  models  are  small  for  step  changes  in  the 
four  considered  inputs.  The  influence  of  Tf  sl  and  TfS3  is  not  shown 
but  has  also  been  evaluated,  and  it  has  been  seen  that  these  tem¬ 
peratures  cause  more  linear  responses  than  Tf,S2  •  Therefore,  from 
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Fig.  2.  Molar  flow  of  hydrogen  over  time  under  incremental  inputs  in  Fc2H50H,in. 
FH2o,in.  Fg  in  and  Tf'S2- 


Fig.  3.  Molar  flow  of  carbon  monoxide  over  time  under  incremental  inputs  in 

Fc2H50H,in.  Fh2 o,in >  Fgiin  and  rf,S2. 

this  comparison  of  the  linear  and  the  non-linear  models,  it  can  be 
concluded  that,  although  it  will  be  important  to  take  care  of  the 
influence  of  Tf  S2,  the  linear  model  can  be  considered  a  valid  con¬ 
trol  analysis  tool.  The  offsets  between  the  linear  model  and  the 
non-linear  model  profiles,  are  indicated  numerically  in  Table  1. 

2.6.  Time  evolution  of  the  output  variables  due  to  input  changes 

The  time  evolution  of  the  reforming  process  variables  is  crit¬ 
ical  because  the  fuel  processor  needs  to  regulate  the  amount  of 
hydrogen  provided  to  the  fuel  cell  stack  (anode)  to  avoid  starva¬ 
tion  or  waste  of  hydrogen  [10].  In  Figs.  4-6,  the  resulting  Fh2  and 
Fco  profiles  at  the  output  of  stages  1 ,  2  and  3  are  plotted.  We  have 


Table  1 

Gap  between  linearized  systems  at  operating  point. 


Linearization  points 

Gap 

OPi 

1.1  x  10“5 

OP2 

2.0  x  10“5 

OP3 

2.0  x  10-6 

OP4 

4.1  x  10-4 

is 

s 


o 

£ 


Si 

& 


Fig.  4.  Disturbances  of  the  system,  based  on  OPn  nominal  conditions  (see  Table 


excited  the  system  by  applying  the  disturbances  represented  in 
Fig.  4. 

Fig.  5  shows  hydrogen  flows  at  the  end  of  each  one  of  the  reform¬ 
ing  process  stages  when  disturbing  the  system  as  indicated  in  Fig.  4. 


5. 


Entrance  Temperature  [K]  £  ,  Fm:  Stage  2  and  3  [mol  min '] 


Table  2 

Input  values  at  nominal  operating  point. 
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Pc2H50H,in  [xlO  3  mol/s] 

pH20,in  [xlO-3  mol/s] 

TS, in  [K] 

rf,Si  [K] 

TfS2  [K] 

T,S3  [K] 

OPn 

1.34 

8.21 

648 

648 

673 

633 

High  magnitude  instantaneous  peaks  occur  because  the  increase  of 
the  volume  flow  at  the  entrance  of  the  reactor  results  in  an  imme¬ 
diate  increase  in  the  flow  at  the  reactor  exit.  For  Fc2h5oh  and  Tgas 
disturbances,  this  peak  is  an  inverse  response.  After  these  peaks, 
with  a  slower  dynamics,  the  newFH2  steady-state  value  is  achieved. 
As  expected,  the  time  required  to  achieve  the  steady-state  follow¬ 
ing  variations  in  FC2h5oh>  Fgas  and  FH2o  is  less  than  that  following 
disturbances  in  Tf  sl,  Ff  S2  and  Tf  s 3.  It  is  also  relevant  to  notice  that 
for  Tf  sl  and  TfS2,  the  time  constant  for  positive  changes  is  smaller 
than  the  time  constant  for  negative  changes.  All  these  observa¬ 
tions  will  be  important  at  the  controllers  design  stage.  In  Fig.  6  the 
behaviour  of  Fc o  is  shown,  which  is  similar  to  that  reported  above 
for  Fh2.  It  is  outstanding  the  influence  of  Tf  S2  at  the  exit  of  stage  2. 

3.  Steady-state  sensitivity  analysis 

In  accordance  to  Figs.  2  and  3,  from  a  steady-state  sensitivity 
analysis  it  is  seen  that  an  increase  of  10%  in  FC2h5oh  with  respect  to 
its  nominal  value  originates  an  increase  in  Fh2  and  Fc o  at  the  outlet 
of  stage  3  of  8%  and  23%,  respectively.  Therefore,  it  is  not  appropri¬ 
ate  to  increase  FC2h5oh  when  more  FH2  is  required  because  larger 
quantities  of  Fc o  will  be  produced  as  well.  On  the  other  hand,  an 
increase  of  1 0%  in  Fh2o  with  respect  to  its  nominal  value  produces  a 
decrease  in  Fh2  and  Fc o  of  3%  and  20%,  respectively.  Therefore,  it  can 
be  said  that  the  sensitivity  of  Fc o  facing  FC2h5oh  and  FH2o  changes 
is  higher  than  the  sensitivity  of  Fh2. 

Concerning  the  effect  of  temperature,  when  the  furnace  tem¬ 
perature  in  zone  1  is  increased  in  10%,  an  increase  in  Fh2  and  Fc o 
of  6%  and  1 7%,  respectively,  takes  place.  When  the  furnace  temper¬ 
ature  in  zone  3  is  increased  in  10%,  an  increase  in  Fc o  of  37%  and 
a  reduction  in  Fh2  of  1.5%  are  observed.  This  is  due  to  the  water 
gas  shift  equilibrium,  which  yields  CO  at  the  expense  of  H2  at  high 
temperatures.  This  input  variable  is  then  very  important  as  control 
variable  because  of  the  sensitivity  of  the  system  to  produce  excess 
CO  when  this  temperature  is  changed. 

Because  of  the  special  behaviour  caused  by  changes  in  the  fur¬ 
nace  temperature  of  stage  2  found  in  the  previous  analysis,  the 
sensitivity  of  the  system  in  front  of  this  variable  is  carefully  con¬ 
sidered.  With  this  aim  we  have  plotted  the  steady  state  values  of 
the  outputs  flowrates  (Fh2  and  Fc 0)  following  changes  in  ?f,S2  with 
respect  to  the  nominal  operating  point  (OPn  in  Table  2). 

In  Figs.  7  and  8  the  hydrogen  yield,  r]H 2,  hydrogen  molar  flow 
rate,  Fh2,  ethanol  conversion,  Xc2h5oh>  anc*  acetaldehyde  conver¬ 
sion,  Xc2h4o>  are  plotted  at  different  7f,S2  values.  From  Fig.  7  it  is 
seen  that  an  increase  of  the  furnace  temperature  in  stage  2  origi¬ 
nates  a  slight  increase  in  Fh2,  whereas  rjH 2  increases  with  a  higher 
rate.  According  to  the  reaction  scheme,  acetaldehyde  conversion 
(*c2h4o)  increases  linearly  with  respect  to  temperature  in  stage  2 
while  the  ethanol  conversion  (xc2h5oh)  is  kept  constant  (Fig.  8). 
Fig.  9  shows  FCo  at  the  outlet  of  stages  2  and  3.  It  can  be  seen  that 
the  influence  of  Tf>s 2  over  the  CO  production  is  very  high. 

4.  Controllability  study 

In  this  section,  a  controllability  analysis  is  performed  using  the 
Relative  Gain  Array  (RGA),  Condition  Number  (CN)  and  Morari 
Resiliency  Index  (MRI). 

One  of  the  most  common  approaches  to  control  a  multiple-input 
multiple-output  (MIMO)  system  is  to  use  a  diagonal  controller, 
which  is  often  referred  to  as  a  decentralized  controller.  The  decen- 


Fig.  7.  Molar  flow  (solid  lines)  and  yield  of  hydrogen  (dotted  lines)  at  different  stage 
2  furnace  temperature  values. 

tralized  control  works  well  if  the  system  is  close  to  diagonal,  which 
means  that  the  plant  can  be  considered  as  a  collection  of  individual 
single-input  single-output  (SISO)  subsystems  that  do  not  interact 
and  can  be  considered  independently.  If  an  off-diagonal  element 
is  large,  then  the  performance  of  the  decentralized  controller  may 
be  poor  [2].  The  controllability  analysis  based  on  RGA,  CN  and  MRI 
permits  to  compare  and  select  the  control  variables  that  minimise 
the  interactions  between  crossed  inputs  and  outputs. 

4. 1 .  Relative  Gain  Array  (RGA ) 

The  Relative  Gain  Array  is  an  analytical  tool  used  to  determine 
the  optimal  control  structure  of  a  MIMO  system.  The  RGA  is  a  nor¬ 
malized  form  of  the  gain  matrix  that  describes  the  interactions 
between  inputs  and  outputs.  Through  the  RGA,  the  process  interac¬ 
tion  of  open-loop  and  closed-loop  control  systems  is  measured  for 
all  possible  input-output  variable  pairings.  A  ratio  between  open- 
loop  gains  to  closed-loop  gains  are  determined  and  the  results  are 


Fig.  8.  Ethanol  conversion  (solid  lines)  and  acetaldehyde  conversion  (dotted  lines) 
at  different  stage  2  furnace  temperature  values. 
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Table  3 

Output  values  at  nominal  operating  point. 


Fh2 ,out  [xlO-3  mol/s] 

Fco.out  [xlO-3  mol/s] 

nH2  [%] 

*C2H5OH  [%] 

*c2h4o  [%] 

yco  [%] 

OPn 

6.39 

1.34 

79.7 

93.4 

84.5 

0.83 

x  ICf4 


Fig.  9.  Molar  flow  of  CO  as  function  of  the  incremental  input  Arf  S3  in  the  outlet  of 
stages  2  and  3. 

displayed  in  the  RGA  matrix.  From  the  RGA  analysis,  different  rules 
can  be  derived  for  the  selection  of  the  appropriate  control  vari¬ 
ables  and  the  best  pairings  between  the  selected  inputs  and  the 
controlled  outputs  [2].  Basically,  RGA  matrixes  close  to  the  identity 
matrix  are  preferred  and  control  structures  with  high  RGA  elements 
should  be  avoided. 

To  perform  the  analysis  of  controllability  of  the  reforming  pro¬ 
cess,  all  possible  combinations  of  two  inputs  among  the  six  possible 
inputs  have  been  considered  (Tables  3  and  4).  The  analysis  is  firstly 
done  at  zero  frequency.  At  this  frequency,  pairs  1,  4,  8,  11  and  13 
are  the  best  because  the  values  of  RGA(  1 ,1 )  are  close  to  one.  In  con¬ 
trast,  pairs  6,  7  and  9  can  be  discarded  because  they  have  negative 
elements  in  the  RGA  at  steady  state.  We  also  eliminated  pairs  2,  3, 
5,  10,  12  and  14  because  with  their  high  values  in  RGA(1,1),  they 
would  be  very  sensitive  to  input  uncertainties. 

4.2.  Condition  number  (CN) 

This  index  is  the  ratio  between  the  maximum  and  minimum 
singular  values  of  the  gain  matrix.  High  CN  values  indicate  that  it 


Table  4 

Controllability  index  RGA,  CN  and  MRI  in  w  =  0  (rad/s). 


PAIR 

RGA(O) 

CN(0) 

MRI(O) 

Pair  1:  Fc2h5oh-Fh2o 

1.40 

7.46 

4.2  x  10-1 

Pair  2:  Fc2H50H-Fgjin 

604 

4.47  x  10“3 

1.2  xl0“3 

Pair  3:  Fc2H5OH-Tf;si 

6.89 

86.40 

1.1  x  10-1 

Pair  4:  Fc2H50H-Ff,s2 

9.5  xlO”1 

5.36 

7.5  x  10-1 

Pair  5:  Fc2H50H~Tfts3 

11.50 

92.11 

2.65xl0-2 

Pair  6:  FH2o-Fgim 

-4x  10-1 

9.66 

0.53 

Pair  7:  FH2o-Ff,si 

-5.1  x  IQ”1 

120.65 

7.8  x  10-2 

Pair  8:  Fu2o~Tfs2 

00 

k 

X 

o 

15.95 

2.56  x  10-1 

Pair  9:  FH2o-Ff,s3 

-0.48 

15.95 

2.56 

Pair  10:  Fg,in-Ff,si 

6.95 

49.59 

2.12  x  10-1 

Pair  11:  Fg,in-Ff,S2 

9.5  xlO”1 

5.16 

1.1 

Pair  12:  Tgtin-T^S3 

11.79 

173.91 

2.75  x  10-2 

Pair  13:  Ff,si-Ffjs2 

9.4  xlO”1 

5.92 

1.33 

Pair  14:  Ff,si-FfS3 

11.79 

173.91 

2.75  x  10-2 

Pair  15:  Tfs2-TfiS3 

5.66  x  10“2 

10.59 

3.34  x  10-1 

will  be  more  difficult  to  control  the  process  because  of  the  sensitiv¬ 
ity  to  uncertainties.  For  this  reason  we  should  select  a  set  of  inputs 
and  outputs  resulting  in  a  system  with  small  CN.  Looking  at  CN 
values  compiled  in  Table  4,  pairs  1,  4,  8, 11  and  13  have  the  lower 
(preferred)  values. 

4.3.  Morari  Resiliency  Index  (MRI) 

The  MRI  (Morari  Resilience  Index)  indicates  whether  a  set  of 
controlled  variables  and  manipulated  variables  provides  a  sim¬ 
ple  control,  giving  a  measure  of  the  inherent  controllability  of 
the  process:  large  values  of  MRI  indicate  that  the  process  is  more 
controllable.  The  pairs  1,  4,  8,  11  and  13  have  also  the  higher 
MRI. 

Having  analysed  the  three  controllability  indexes  at  steady- 
state,  it  can  be  concluded  that  the  control  structures  consisting  in 
pairs  1 , 4, 8, 1 1  and  1 3  are  better  than  the  other  ones.  For  this  reason, 
in  the  next  section,  only  these  structures  will  be  considered. 

5.  Selection  of  the  control  structures 

In  this  section,  the  steady  state  controllability  analysis  is  com¬ 
pleted  with  the  analysis  at  different  frequencies  [2],  and  the  best 
control  structures  are  identified.  Only  pairs  1,  4,  8,  11  and  13  are 
taken  into  account.  Pair  1  corresponds  to  of  the  manipulation  of 
Fc2h5oh  for  the  control  of  Fh2,  and  the  manipulation  of  FH2o  for  the 
control  of  Fc o-  On  the  other  hand,  for  the  pairs  4,  8,  11  and  13  if 
we  need  to  control  the  Fh2  we  should  act  on  FC2h5oh»  FH2o*  Tgin  and 
7>,si,  respectively.  To  control  the  Fc o  we  must  act  on  the  7f,S2- 

5.1.  Controllability  Sequential  analysis 

The  frequential  analysis  is  done  considering  the  frequency  range 
from  1 0-3  to  1 02  rad  s-1 .  In  Figs.  10-12  the  results  obtained  for  the 
different  controllability  indexes  are  shown.  In  the  three  figures,  the 
order  of  the  different  pairs  is  maintained  during  a  wide  frequency 
range.  The  peaks  observed  in  the  three  figures  are  due  to  numerical 
problems  and  occur  at  different  frequencies  if  the  system  model  is 
trunked  with  a  different  number  of  states.  Pairs  4,  11  and  13  are 


V.M.  Garda  et  al.  /  Journal  of  Power  Sources  196  (201 1)4411  -4417 


4417 


Fig.  11.  Morari  Resiliency  Index  (MRI)  of  different  pairs. 


Frecuency  [rad/s] 


Fig.  12.  Condition  number  (CN)  of  different  pairs. 

the  best  in  accordance  to  the  three  controllability  indexes.  Their 
RGA(1,1 )  and  MRI  (close  to  one)  and  CN  (lower  than  10)  are  accept¬ 
able  values  for  scaled  systems.  In  the  three  cases,  to  control  the  Fc o 
the  preferred  manipulated  variable  is  Tf  S2. 


6.  Conclusions 

This  work  focuses  on  the  design  of  controllers  for  an  ethanol 
steam  reformer,  which  is  a  MIMO  system  with  six  inputs  and 
two  controlled  outputs.  The  main  tasks  of  the  work  are  the  char¬ 
acterisation  of  the  dynamic  response  and  the  selection  of  the 
preferred  control  structures.  Both  tasks  are  preliminary  studies 
for  the  design  of  controllers.  The  dynamic  response  is  analysed 
through  a  non-linear  model.  Inverse  responses  and  non-linearities 
are  observed.  The  controllability  analysis  is  based  on  a  linear  model. 
In  accordance  to  RGA,  MRI  and  CN  controllability  indexes,  three 
different  control  structures  are  selected  as  the  most  promising 
ones  for  the  control  of  Fh2  and  Fc o-  They  are  the  pairs  con¬ 
sisting  in  the  manipulation  of  Fc2h5oh  and  TfS2,  Tg in  and  TfS2, 
and  Tf.si  and  Tf,S2-  All  three  pairs  include  the  Tf,S2-  These  control 
structures  are  a  priori  appropriate  for  2  x  2  MIMO  decentralised 
control.  However,  due  to  the  non-linearities  of  the  system,  spe¬ 
cially  apparent  in  the  rfS2  response,  the  performance  of  the 
controllers  cannot  be  confirmed  until  its  validation  in  non-linear 
models. 
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